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Abstract 

A recently developed procedure for a partial wave decomposition of a three-nucleon force is 
applied to the tt — vr , tt — p and p — p components of the Tucson-Melbourne three-nucleon po- 
tential. The resulting matrix elements for the vr — vr and tt — p components are compared with 
the values obtained using the standard approach to the partial wave decomposition, in which the 
TT — p expressions for the matrix elements are also derived and presented. Several numerical tests 
and results for the triton binding energy and the correlation function prove the reliability and 
efficiency of the new method. 



PACS numbers: 21.45.-v, 21.30.-x, 21.45.Ff 



I. INTRODUCTION 



The Tucson- Melbourne (TM) three nucleon force (3NF) is an important model of the 
three-nucleon (3N) interaction. It consists of three parts stemming from exchanges of tt — tt , 
TT — p and p — p mesons. The main ingredient of the TM force, the meson-nucleon scattering 
amplitude with the off-shell mesons, was derived using the current algebra techniques. This 
was done in [l[ and improved in [3] for the tt — tt part. The n — p and p — p contributions 
were derived in j2l, 0] . In jsf the structure of the tt — tt part of the TM 3NF was revisited to 
achieve a consistency with the chiral symmetry and the modified force is known as the TM' 
model. 

The effects of all terms on the triton binding energy were studied in [6|. It turned 
out that the n — p force acts repulsively for the '^H contrarily to the n — n interaction 
and combining them leads to the binding energy close to the experimental value. The 
p — p force has only a small influence on the triton binding energy. A similar behaviour 
was observed for scattering observables in the three-nucleon system the largest effects 
came from the dominant tt — vr part and the influence of the n — p part was smaller and in 
the opposite direction. The p — p contribution proved to be much smaller and practically 
negligible. However, the results of Refs. j^, 0] were based on partial waves restricted to 
the total angular momenta in the two- nucleon subsystem j < 2. Thus conclusions of 0] 
are valid only in a low energy domain of the three-nucleon continuum. For higher energies, 
where more partial waves are required to achieve convergence, only the tt — tt part of the TM 
force was used, (see e.g. js], 0]). While the inclusion of this main component of the TM 3NF 
improves the description of many scattering observables, some serious discrepancies with 
data remain and they become larger at higher energies. One of the possible explanations 
for this disagreement is a lack of shorter-range parts of the 3NF in those calculations, what 
calls for a reliable and fast method to obtain matrix elements for all components of the TM 
force in higher partial waves. 

Recently we have proposed a novel, automatized way to perform a partial wave decom- 
position of any two- and three-nucleon potential jlo| . This approach makes use of a software 
for symbolic calculations to generate the part of the code which is specific for a considered 
force model. More precisely, in this way we calculate exactly the isospin and spin-momentum 
parts of the nuclear interactions and generate a corresponding FORTRAN (or C) code. That 
momentum dependent output forms an integrand for further five- dimensional numerical in- 
tegrations. 

In this paper we present results of applying that new scheme to the original TM 3N force. 
They confirm the feasibility and efficiency of our method and its numerical implementation. 
The existence of such a reliable procedure is especially important in view of available and 
forthcoming results from the chiral perturbation theory (xPT) for 3N forces at higher 
orders of the chiral expansion. A big number of different momentum-spin-isospin structures 
contained in those interactions requires a safe and automatized method to perform partial 
wave decompositions, which is guaranteed by our method. Furthermore our scheme avoids 
the application of partial wave decomposed permutation operators when dealing with prod- 
ucts of 3NF's and permutation operators as they are often required, e.g., in 3N Faddeev 
equations. Such an application is numerically demanding because it requires a huge number 
of partial waves. Thus, again an efficient, fast and precise method is needed. 

Our novel scheme of an automatized partial wave decomposition (aPWD) is described in 
Sec. [TTl Results and additional tests for our numerical realization are presented in Sec. IIIII 
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and conclusions are given in Section IIVI The standard PWD of the tt — p component of the 
TM 3NF is given in the Appendix. 



II. AUTOMATIZED PARTIAL WAVE DECOMPOSITION 

The 3NF, V123, is an indispensable ingredient in a theoretical description of the few-body 
systems. It can be always written as a sum of three terms 

V,2s = V^'^ + V^(2) + V^^^ , (2.1) 

where each jg symmetrical under the exchange of nucleons j and k {i,j,k = 1,2,3, 
i j k). Such a splitting in the case of the vr — tt exchange TM 3NF corresponds to the 
possible choices of the nucleon undergoing off-shell vrN scattering. 

The 3NF typically enters the dynamical equations via its part V^^\ In the case of the 
three-nucleon bound state, the Faddeev component ifj fulfils the following equation 13 



^ = GotP^ + (1 + Got)G'o V^(^) (1 + P)V' , (2.2) 

where Go is the free 3N propagator and t is the two-body t-operator generated from a 
given nucleon-nucleon (NN) potential through the Lippmann-Schwinger equation. The per- 
mutation operator P = P12P23 + -P13-P23 is given in terms of the transpositions Pij, which 
interchange particles i and j. The full bound state wave function is then obtained as 
^ = {1 + P)tfj. 

Transition amplitudes for the elastic nucleon-deuteron scattering, U, and for the breakup 
reaction, Uq, are given as [l3| 



U = PGo ^$ + PT + -1/^(1 + P)$ + r(^)(l + P)G'oT , 
Uo = {1 + P)T , (2.3) 

where the auxiliary state T fulfils the 3N Faddeev equation 

T = tP$ + (1 + tGQ)V^^\l + P)^ + tPGoT + (1 + tGo)V^^\l + P)Gor , (2.4) 

with $ being the initial state composed of the deuteron wave function and a momentum 
eigenstate of the projectile nucleon. 



states I g, a) in the j J-couphng [ij, ll5| 



Equations (12. 2p and (12.41) are solved ll3lJl4| in the momentum space using 3N partial- wave 



I p,g,a) ^1 pq{ls)j{\]^)I{jI)JMj) \ (t^)TMT) , (2.5) 

where p and q are magnitudes of the standard Jacobi momenta and a denotes a set of discrete 
quantum numbers arising in the following way: the spin s of the subsystem composed from 
nucleons 2 and 3 is coupled with their orbital angular momentum I to the total angular 
momentum j. The spin | of the spectator particle 1 couples with its relative orbital angular 
momentum A to the total angular momentum of nucleon 1, J. Finally, j and / are coupled 
to the total 3N angular momentum J with the projection Mj. For the isospin part, the 
total isospin t of the (23) subsystem is coupled with the isospin | of the spectator nucleon 
to the total 3N isospin T with the projection Mp. 
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Any three-nucleon force enters Eqs. f|2.2p - fl2.4l) in the form of V^^\l + P). Therefore a 
partial wave decomposition of V^^^ as well as V^^^P has to be performed. The standard 
approach to perform a partial wave decomposition of V^^'' is very tedious, even with 



improvements suggested in [17[, since each momentum-spin- isospin structure, which occurs 
in a 3NF, has to be treated separately. In the case when a 3NF consists of a big number 
of such structures, like chiral 3NF's at higher orders of the chiral expansion, the traditional 
approach to a partial wave decomposition is very inefficient and extremely time consuming. 
In addition, the application of the permutation operator, when calculating V^^^P, causes an 
additional numerical problem, which originates from a slow convergence of the V^^'^P matrix 
elements with respect to the number of intermediate states \a" >: 

{p,q,a\V^'^P\p',q',a')= (2.6) 
= I dp"p"' I dq"q"' q, a \ V^''> \ p" , q" , a") (/, q" , a" \ P \ p' , q' , a') . 

a" 

In order to calculate precisely these matrix elements, a big number of intermediate states 
I a") is required, and, thus, one is forced to calculate matrix elements of V^^^ operator for a 
much bigger set of a" states than actually needed in order to get converged solutions of the 
Faddeev equations. 

In our new approach, called in the following automatized partial wave decomposition 
(aPWD), to get matrix elements of V^^^ and V^^^P, that drawback is removed because 
matrix elements of V^^) and V^^^P are calculated directly. 

The starting point of our method is the observation, that any 3N interaction and thus 
also its V^*^^-* component in momentum space can be written as a sum of terms in the form 

^^^^ = fiol, 02, q3)0spin{qi, 02, 03, 0^1, ^2, (^3)Oisospin{'n, fg, fg) (2.7) 

where Ospin and Oisospin are operators acting on spin and isospin degrees of freedom, respec- 
tively, which are built from the spin (cij) and isospin (r^) operators of individual nucleons. 

Scalar factors f{ql, 02,03) and spin operators Ospmiqi, (h, Os, (^1, '^2, o's) depend on momen- 
tum transfers ql to nucleon i which are expressed in terms of the initial and final Jacobi 
momenta p, q and p', q', respectively, as: 

qi = q' -q, <i2 = {p' -p)- ^(g' - g) , 93 = -{p' -p) -^{q' -q) = -{qi + & ■ (2.8) 

For example, in the vr — tt part of the TM 3N force one meets the following spin-isospin 
structures: 

Ospiniqi, q2, qs, o\, (^2, (Js) = (0^2 ■ & {(Ts ■ (B) , (<T2 ■ g"i) ((fg ■ ^3) (g^ ■ gi) , 

{(^2 ■ & {(Ts ■ qa) {{q2? + (qsY) , cfi ■ q2 x ^3 (2.9) 

Note, that not all combinations of Ogpin and Oisospin actually appear in the above example. 

In the first step of aPWD we calculate 3NF matrix elements using partial wave states 
I P, q, P) pJi] in the so-called L6'-coupling 

I P, q, P)i =1 pq{l\)L{s]^)S{LS)JMj), \ {t\)TMT)i , (2.10) 
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where the relative orbital angular momentum / (within the pair (23)) and A (between the 
pair (23) and nucleon 1) are coupled to the total orbital angular momentum L. In the spin 
space, the spin of the (23) pair is coupled with the spin ^ of the nucleon 1 to the total spin 
S. Finally L and S are coupled to the total 3N angular momentum J with the projection 
Mj. The index 1 emphasizes, that the spectator particle is nucleon 1. /? describes the set 
of discrete quantum numbers discussed above. The isospin state is the same as in the basis 
state \p, q, a). 

In this basis, it is easy to decouple the isospin and spin parts from the momentum part, 
what leads to the following form of a 3NF matrix element 

{p'q'{l'X')L'{s'^)S'{L'S')JMj I {{t'^)T'mT' \ V'^'^ \ pq{lX)L{s^)S{LS)JMj) \ {t^)TMT) 
= fdp'f dq'j dp! dqY,C{L\S\J-mL',Mj-mL',Mj)Y,C{L,S,J-mL,Mj-mL,Mj) 

mj^i mi, 



X f{p\q',p,q){{t'\)T'MT \ O.^ospin \ {t\)TMT), (2.11) 
where 

I 

yi^r'^iP^ ^) = X] '^(^' ^' '^l^ - mi, mi) Yi^rniiP) Yx,mL-m,iq) ■ (2.12) 

■mi=-l 

with the standard Clebsch-Gordan coefficients and the spherical harmonics. For abbreviation 
we skip in (12. lip and in the following the spin (fi and the isospin fi operators in the arguments 

of Ospin and Oisospin operators. 

The matrix element in the spin space appearing in (12. lip . 

{{s'^)S'Mj — niL' I Ospinip', q',P, q) \ {,s\)S Mj — rriL), depends on the momenta qi and s pin 



quantum numbers. Using a software for symbolic calculations (such as Mathematica® 18 
in our case) it is very easy to calculate this matrix element for all combinations of spin 
quantum numbers as a function of the momentum vectors q^. To this aim we use the 
Kronecker product built in Mathematica, which allows us to express the spin matrix 
element in terms of simple matrix operations. This is even more straightforward in the 
case of the isospin matrix element, which does not depend on any additional parameters. 
Another advantage of using software for symbolic calculations is the possibility to generate 
a Fortran (or C) code in an automatized way. This eliminates possible errors which can be 
introduced during programing of very lengthy formulas for the spin matrix element. The 
calculation of the 3NF matrix elements requires finally an eight-dimensional integration 
shown in (12. lip . In a typical case the total isospin and its projection is conserved. We also 
assume that the considered 3N force is rotationally invariant. Then the matrix elements 
in (12. lip vanish unless J = J' and Mj = Mj', and, additionally, do not depend on Mj. 
Thus we can calculate 

G{1', A', L', s', S", t', I, A, L, s, S, J, t, T, Mt) = 

■^1 1 1 



J2 {{t'-)T,MT I {pW{iy)L'{s'-)S\L'S')JMj I I pq{lX)L{s-)S{LS)JMj) 

Mj=-J 



X I (4)T,Mt), (2.13) 
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which is equal to the original matrix element of V^^^ given in Eq. fl2.1ip . The integrand in 
G{1', A', L', s', S', t', I, A, L, s, S, J, t, T, Mt): 



X 



CiL, S, J; m^, Mj-m^, Mj)y;';r'' {p' , q') y^^^ip, q) 



X {{s'l)S'Mj - rriL' \ Ospinip', q'',p, q) \ isl)S Mj - m^) 



X f{p',q'',p,q){{t'^)T'MT \ 6..o.p^n I (4)™t) (2.14) 



is a scalar and thus does not depend on all directions of the Jacobi momenta p^]. Therefore 
we are free to choose for example p along z-axis {p = (0, 0,p)) and (pg = and are left with 
five-fold integrations only 

G{1', X', L', s', S', t', I, A, L, s, S, J, t, T, Mt) 
= 8n' I dp' I dq'l dcosiOg)-^ ^ Y.'^^^' ^ '^'^'^^'^^J ' "^^'^^j) 

J J J + M.j=-Jm^, 

Y cm S, J; m^, Mj - tul, Mj)3^; J''"^' (p', q') 3^i"^^(5, q = (sin(^,), 0, cos(^,))) 



X 

rriL 



X {{s'^)S'Mj-mL' I dspin{p',q'',P= (0, 0,p), g = g(sin(^,), 0, cos(^,))) | {s^)SMj-mL) 
X f{p'A'J= (0,0,p),g = g(sin(^,),0,cos(^,)))((t'i)TMT | ^,sosp^n I (2.15) 

The reduction of the number of integrations for a simple example of 3NF is numerically 
exemplified in Ref. 10 1. 

The remaining summations over m^/ , and Mj and 5-fold integrations can be performed 
for a small number of (p,q,p',q') combinations even on a personal computer. However, 
a large number of five-dimensional integrations, as they are needed to obtain all matrix 
elements needed for the solution of the 3N Faddeev equations, has to be carried out on a 
powerful parallel computer. Once the matrix elements {p' , q' , /3' \ V^^^ \ p, q, (3) are calculated, 
recoupling to the ^'/-representation, {p' , g', a' \ y*-^^ | p, g, a), can be easily performed [l5| 

[ I s J 

(p', g', a' I V^^^ \p^q,a) = Y, V(2j + 1) (2 J + 1) (2L + 1) {2S + 1) I A \ I 

13,13' \^ L S J 

( I' s' f ] 

X V(2j' + l)(2J' + l)(2L' + l) {2S' + 1){ X' I r Wp',g',/3'|l^W|p,g,/3).(2.16) 

[l' S' J ) 

Now let us turn to the ^''^•'(l -|- P) operator and discuss its V^^^ P12P23 matrix element 
i(p',g',/3' I y«Pi2P23 I P,g,/3)i = [ dp' [ dq' [ dp [ d^ i(p',g',/3' | p'q'){p'q' \ y« 



X Pi2P23\pq){pq\p,q,l3)i. (2.17) 
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Since 



P12P23 I = I -\p+ \l -P- l^iP^Pr I {s\)SMs), 

^ pi^ospinpisospin | (^t-)TMT)l (2.18) 

PttPr I {sl)SMsh = I (4)'5M5)2 = 

= J](-)^V^|| '^'} I {s"l)SMs), (2.19) 



= \ I (i"^)™T)i , (2.20) 

where p^P^'^(^p^^°^P'''^^ jg the part of operator acting in the spin (isospin) space, one gets 
i{p',q',/3' \V^'^PuP23\p,q,l3)i^ J dp'J dq'J dpi dqJ2C{L',S',J;mL',Mj-mL',Mj) 



X i{{s'^)S'Mj - mv I dspin{^, q', ~p+ ^q, -p- ^q) \ {s"^)S Mj - mL)i 
X /(p', q', -\p+ \q, -P- \q)i{{t']^)T'MT \ A.o.pm I (t"^)™T)i . (2.21) 
Similarly, for ^'^^^^13^23 one gets 

I F(')Pi3P23 I dp'l dq'l dp! dqJ2CiL',S',J;mL',Mj-mL',Mj) 

X ^C{L,S,J-mL,Mj-mL,Mj) 

rriL 

X y;;(r-{j/,q')ytr{p,q)Y:i-r^s^^^ \ i '1} 



2 " ■> t" 2 

X i((s'^)5'Mj - niL' I api„(p'', g^', ~p- \l p- \q) \ {s"]^)S Mj - mL)i 

X /(p', q\ -\p- \qj- \q),{{t']^)T' Mt \ O.sosmn \ (t"^)™^)!. (2-22) 

That means that the calculation of these two contributions proceeds in the same way as 
calculation of the V^^^ matrix element. Only the arguments of the term Ospin have to 
be changed and additional factors originating from the recoupling of the spin and isospin 
quantum numbers have to be taken into account. As for the V^^^ operator also here the 
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TABLE I: The values of the discrete quantum numbers for selected a states ()2.5p for the total 
angular momentum J = \ and the positive parity 11 = (—1)'"'"'^. 



eight-fold integrations can be reduced to the five-fold ones and recalculation to |p, g, a) 
states can be performed. 

It is important to note that, since our basis states g, a) are antisymmetric with respect 
to the exchange of nucleons 2 and 3, (2.19) and (2.20) yield the same values for the matrix 
elements. This allows one to reduce significantly the size of the codes and the required 
computation time. 



III. RESULTS 

A. The TM 3NF and its tt — vr , vr — p , and p — p components 

Since the aim of this work is not to study the dependence of the matrix elements of the 
TM force on its parameters, in the following we use their values given in Tab. I of Ref. 
a = 1.03^~\6 = -2.62^-^c = 0.91^~^(^ = -0.753/i-3 with ^ = 139.6 MeV. In the 
numerical implementation of f l2.1ip we use the same number of gaussian points for each 
of the five angular domains. It might be more efficient to relax this constraint in future 
applications and to optimize the grids further. Thus our integration method leaves room 
for improvement, even if we will later demonstrate in subsection IIII Gl that it leads to fully 
converged results. 

The TM 3NF matrix elements calculated in the basis (12. 5p are functions of four momen- 
tum magnitudes and two sets of discrete quantum numbers. In Figs. [I]l2l examples of the 
TM force V^-^^ matrix elements are shown together with its tt — tt, tt — p and p — p compo- 
nents in one- dimensional plots. In Fig. [H the matrix elements < p', g', a'|y''^^|p, g, a > for 
p' = q' = q = 0.132 fm~^ and for different channel pairs (a', a) (see Tab. [I]) are shown as a 
function of the momentum p. The same matrix elements but for the momenta p' = 0.711 
fm~^, q' = 0.132 fm~^, and q = 2.84 fm~^ are shown in Fig. [2] again as a function of p. The 
71 — n part dominates in all cases but the n — p part is also important (see Fig. [T)Dj2^-c). 
The p — p part is of less importance for all the considered matrix elements. 



B. The aPWD for V(^)(l-|-P) operator 

As was described in Sec. [Tll aPWD can be 
to the y^^)(l -I- P) operator. Using aPWD for 
the V^^'' operator: the automatized procedure 
reduces the possibility of errors. In the current 



applied not only to the V^^'^ alone but also 
\/(i)(l + P) has the same advantages as for 
can be easy tuned to any kind of 3NF and 
implementation of aPWD the calculation of 
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<p',q',a' I V^^^ I p,q,a > [fm^] 
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FIG. 1: (Color online) The TM 3NF matrix elements < p' = 0.132/m-\g' = 
0.132fm~^,a'\V^^^\p,q = 0.132/m~^,a > as a function of the p momentum for (a', a): a) (1,1), 
b) (1,4), c) (6,3) d) (6,8). The solid (red) curve represents the full TM 3NF and the other curves 
show the contributions coming from the vr — vr (black dotted) , n — p (blue dashed) and p — p (green 
dot-dashed) components. 



V^^\l + P) matrix elements needs about the one and half amount of the computing time 
needed for V^^^ , what is important from the practical point of view. Finally, in the standard 
scheme of PWD, the number of intermediate partial waves used to represent the P operator 
is limited and might be insufficient. In the case of aPWD there is no separate decomposition 
of the permutation operator what corresponds to the inclusion of all three-body intermediate 
waves. In Figs. E] and H] the matrix elements of V*^^-*(H-P) are shown for the same momenta 
and channels as in Figs. [Hand [21 respectively. For the channel combinations (1,1) and (6,3) 
in Fig. [3] and (1,1) and (6,8) in Fig. [U where vr — tt force dominates, the picture is similar 
to the corresponding ones in Figs. [1] and [2l For the remaining channel combinations the 
differences are more visible, for example the inclusion of the permutation operator for the 
TT — vr component for the (6,8) pair in Fig. [3] leads to the change of the sign and strength of 
this force. In that case also the tt — p part becomes bigger after the permutation operator 
is applied. Also for the (6,3) case in Figs. [3] and [H the action of the permutation operator 
changes the strength of the matrix element and increases the momentum range, where both 
TT — TT and 71 — p components play a significant role. For all the here presented cases the 
p — p force is much smaller than the remaining interactions. 
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FIG. 2: (Color online) The same as in Fig. [T] but for momenta values: "p' = 0.711/m ^,q' 
0.132/m-^g = 2.842 fm'^ 



The aPWD method allows us to study the role played by different isospin structures 
entering the TM force. An example is given in Fig. [5] where, for the tt — p force, the 
contribution from the so-called " Kroll-Ruderman" and two "A" terms 0] (see also Appendix 
lA ip are shown. For the presented matrix elements (< p' = 0.132/m~^, q' = 0.132/m~^, a' = 



l\v}]}p{l + P)\p,q = 0.132/m~^,a = 1 >) the "Kroll-Ruderman" term dominates for small 
momenta p, while the two "A" terms are bigger for p > 2 fm~^. However, they have 
opposite signs, so their combined effect is weak and leads to a reduction of the strength of 
the dominant "Kroll-Ruderman" term. 



C. The comparison of the standard and automatized PWD schemes for vr — vr and 
TT — p forces. 

For the vr — vr force the partial wave decomposition has been presented in jsj and in an 
alternative way in [l^ . The comparison of results obtained by aPWD and the ones obtained 
in Ref. [l6| is presented in Fig. O Again the channel pairs and momenta are chosen as in 
Fig. [TJ A very good agreement between the both methods is clearly seen. 

In Appendix A we present expressions for the partial wave decomposition of the tt — 
p force. This decomposition is in the spirit of the decomposition of the tt— vr interaction given 
in Ref. [l6|. In Fig. [7] we compare results obtained in the aPWD scheme with those based on 
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FIG. 3: (Color online) The same as in Fig. [T]but for the + P) operator. 



PWD given in Appendix A. Because of the internal construction of the PWD from Appendix 
A, we compare matrix elements of V'^^^ PxzPiz instead of V'^^^ . The matrix elements of the 
standard PWD are obtained using partial waves up to jmax = 5 in intermediate states. For 
this truncation, the matrix elements considered here are converged (see Sec. IIII El) . Again, 
for all given examples, the agreement between both methods is excellent. 

Though the CPU time is smaller for the scheme presented in Appendix A, the long time 
which is needed for the derivation of the partial wave decomposition of complicated spin- 
momentum structures and its programming in the standard way is incomparable with the 
relatively short time demanded by aPWD. Another advantage of aPWD lies in its flexibility 
which allows one to use it easily for different operators. In the case of the standard PWD 
each spin-momentum structure has to be treated separately. 



D. The equality of y(^)Pi2i^23 and y(^)Pi3P23 

The equality of V^''^^Pi3P23 and V^''^^Pi2-P23 matrix elements between the states antisym- 
metrized in the (23) subsystem forms another nontrivial test of numerics. To check this, 
we compare some matrix elements for V^^^''Pi2-P23 obtained via Eg. (12.211) with the corre- 
sponding ones for V^*^^^Pi3P23 from Eq. (l2.22p . Results are displayed in FigJH] again for four 
combinations of channel pairs and selected values of p', g' and q momenta (the same as in 
Fig. [T]). The numerical confirmation of the equality of the V^*^^^Pi3-P23 and y'^^^P\2P2z matrix 
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FIG. 4: (Color online) The same as in Fig. [2] but for the + P) operator. 



elements is clear. They differ from the y*-^^ elements, as can be seen for some examples in 
Fig. m All three possibilities are shown: for the channel combinations (1,1) and (1,4) V^^"* 
dominates, while V^^'' P12P23 and V^^^Pi3P23 are much smaller. For the (6,3) combination 
each operator gives a similar contribution to V^*-^^(l + P). For the (6,8) choice and momenta 
around 2 fm~^ the contribution from V^*^^) is much smaller than the remaining two. 



E. The convergence of ^(^^(1 + P) matrix elements with respect to the number of 
the intermediate partial waves for the tt — p and the full TM forces. 

The aPWD result for the ^^'^^•'(l + P) operator, which corresponds to the infinite number 
of the intermediate partial waves taken into account during the action of the permutation 
operator, gives the limit to which results of the traditional scheme should converge. This 
convergence is confirmed in Figs. [9] and [10] for the n — p part of the TM and the full TM 
3NF, respectively. The channels and momenta are the same as in Fig. [2l While for the 
channel combination (1,1) already the smallest number of partial waves gives the aPWD 
limit, for the other combinations much more partial waves have to be taken into account. 
For one of the shown here cases (Fig. [TUb ) taking all partial waves up to jmax = 5 is still 
insufficient to achieve the limit of aPWD. Note, however, that the magnitude of this matrix 
element is relatively small. In general, the convergence of the traditional PWD scheme is 
fully confirmed. 
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<p',q',a'l V^ip(l+P)l p,q,a> [fm] 
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FIG. 5: (Color online) The contributions from the different parts of the n — p force for matrix 
elements < p' = 0.132 fm'^g' = 0.132 fm'^.a' = l|V;Lp(l + P)\p,q = 0.132 im-\a = 1 >. The 
black solid line represents the total tt — p TM 3NF while the red dotted, green dashed and blue 
dot-dashed lines represent the " Kroll-Ruderman" , the isospin even A and the isospin odd A terms, 
respectively. 

F. The binding energy and correlation function for ^H. 

As a first application, we would like to calculate in the following the binding energy of ^H, 
some energy expectation values and the correlation function. The obtained binding energies 
and expectation values of the kinetic energy {Hq), the NN potential energy (V/vat) and the 
3N potential energy (V^n) are given in Tab. [ITlfor several realistic NN interactions alone and 
together with the TM force. The TM force was included for all states with subsystem total 
angular momentum j < 2. The inclusion of the TM force leads to a stronger binding of ^H. 
The binding energy changes, after the inclusion of the TM force, by approximately -1.093 
MeV for the CDBonn potential and from -1.122 to -1.334 MeV for Nijmegen potentials. 
These results are in a reasonable agreement with the ones given in Tab. 2 of Ref. for 
the Bonn OBEPQ (-9.596 MeV) and the Nijmegen (-8.689 MeV) potentials. Note, that in 
Ref. slightly different values of the a, b and c parameters were used. In our calculations, 
we include partial waves up to jmax = 5 for the two-body interaction. This is also different 
from Ref. p where only partial waves up to jmax = 2 were included. Of course, for the 
given set of parameters, the binding energies do not accurately reproduce the experimental 
value of -8.482 MeV. Because of the well-known scaling behavior of many N-d scattering 
observables with the triton binding energy (see for example [sO]), it will be necessary to 
finetune the TM model such that the triton binding energy is more accurately reproduced, 
e.g. along the lines of Ref. |1^ 

The inclusion of the TM 3NF leads for all the NN potentials to higher expectation values 
of the kinetic energy and lower expectation values of the NN potential energy (about 3-6 
MeV). The expectation values of the 3N potential energy amounts from 3.3% to 4.8% of 
the expectation values of the NN potential, depending on the particular NN potential. This 
observations are in line with the general expectations for the strength of 3NF's and the more 
compact state of '^H when the binding energy is increased. 
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<p',q',a' I V^^^^ I p,q,a > [fm^] 
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FIG. 6: (Color online) The comparison of the matrix elements of vr — tt force obtained in the 
standard (crosses) and automatized (solid line) PWD. The channel combinations and momentum 
values are the same as in Fig. [H 



The correlation function is defined in the configuration space as l3 

C{r) = ^^J dri"^ I J2Sir-r7j) \ ^) 



(3.1) 



where rjj is the relative distance operator conjugate to the operator of the Jacobi momentum 
p. It is shown in Fig. [11] for the different NN potentials alone and combined with the 
TM 3NF. For the smaller distances shown in Fig. [TTl the probability to find two nucleons 
increases when the TM 3NF is included. At least in part, this can be understood because 
the correlation functions drop more quickly for larger r due to the increased binding energy. 
Note that at short distances, the effect of the 3NF's is much larger than the dependence on 
the NN interaction model. The here presented correlation functions are in good agreement 
with the ones presented in 12|] for the same NN potentials combined with the tt — tt part 
of the TM force. 



G. The quality of the five-dimensional integration 

Finally, we would like to give an example of the stability of aPWD against the number of 
points used in the numerical integrations. In Tab. IIIII the ^'^^^(1 + P) matrix elements are 
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<p',q',a' I V^^ip P^3P23 I p,q,a > [fm'] 
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FIG. 7: (Color online) The comparison of the matrix elements of the tt — p force obtained in the 
standard PWD from Appendix A (crosses) and automatized (solid line) PWD. The channels pairs 
and momenta are the same as in Fig. [2l 



NN potential 


Et [MeV] 


{Ho) [MeV] 


{Vnn) [MeV] 


{Vsn) [MeV] 


CDBonn 


-8.008 


37.620 


-45.609 




Nijmegen I 


-7.738 


40.737 


-48.467 




Nijmegen II 


-7.658 


47.526 


-55.176 




Nijmegen 93 


-7.664 


45.617 


-53.283 




CDBonn + TM 


-9.101 


41.934 


-48.669 


-2.342 


Nijmegen I + TM 


-8.860 


45.523 


-52.277 


-2.098 


Nijmegen II + TM 


-8.992 


54.318 


-61.112 


-2.189 


Nijmegen 93 + TM 


-8.841 


51.173 


-58.092 


-1.925 



TABLE II: The triton binding energies Et and the energy expectation values (Hq), {Vnn) 
(Vsn) for the different NN potentials alone and together with the TM 3NF. 



and 



given for the same channels and momenta as in Fig. [T] for two values of momentum p=0.711 
fm~^ and p=5.959 fm~^. Results were obtained using N=12 or N=15 gaussian points in each 
of the five integrations in (12. lip . The agreement seen in Tab. Illll between both predictions 
clearly demonstrates that the numerical integration is well under control and leads to fully 
converged numbers. 
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FIG. 8: (Color online) The contribution of Oi = V'^^'^PxzP^z and V'^^^Px^P^z operators to the 
total ^'■^•'(1 + P) TM 3NF matrix elements. The channels combinations and momenta are chosen 
as in Fig. [TJ The crosses represent T^^^^i-'i3-P23 matrix elements. The dashed, dotted and solid lines 
represent ^^^^-Pi2-P23 and y*^^)(l + P) matrix elements, respectively. 



(a', a) 


N=12 
p=0.711 fm"^ 


N=15 
p=0.711 fm"^ 


N=12 
p=5.959fm-i 


N=15 
p=5.959fm"i 


(1,1) 
(1,4) 
(6,3) 
(6,8) 


-0.0010139197 
-0.00083291615 
-2.8739479x10"°^ 
-1.212581x10"°^ 


-0.0010139197 
-0.00083291615 
-2.8739479x10"°^ 
-1.21258x10"°*^ 


8.9762335 xl0-°5 
-9.0900123x10"°° 
1.6945111x10"°° 
-3.5797036x10"°^ 


8.9762335 xl0"°5 
-9.0900124x10"°° 
1.6945111x10"°° 
-3.5797045x10"°"^ 



TABLE III: The ^'■^•'(1 + P) matrix elements (in fm^) for channels combinations and momenta 
p', g', q as in Fig. [H depending on the number of gaussian points N used in the five-fold integration 
of Eq. (|2.14|) . The value of momentum p is 0.711 fm"^ (left) and 5.959 fm"^ (right). 



IV. SUMMARY 

We apply an automatized method of partial wave decomposition to the Tucson-Melbourne 
three-nucleon force. The obtained results agree very well with the traditional way of a partial 
wave decomposition for tt — vr and vr — p contributions to the TM 3NF. For the latter one, 
we also give formulas of the partial wave decomposition in the traditional approach. Matrix 
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<p',q',a' I V^^^ (1+P) I p,q,a > [fm^] 
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FIG. 9: (Color online) The convergence of the matrix elements of the n — p part of the TM force: 
V^^p{l + P) with respect to the number of the intermediate partial waves used during the action 
of the permutation operator (see Eq. (j2.7p ). The channel combinations and momenta are the same 
as in Fig. [2l The crosses represent predictions obtained within the aPWD approach, the dotted 
(black), dash-double dotted (red), dash-dotted (green) and solid (black) lines represent results 
obtained with the traditional method described in Appendix A with all the intermediate 3N states 
up to jmax = 2, 3, 4 and 5, respectively. 



elements obtained in the new way are used in the calculations of the triton wave function 
with different underlying nucleon-nucleon potentials. We performed also different numerical 
tests, which confirm the reliability of our method and computer codes. 

Among many advantages of aPWD, we would like to emphasize its generality, efficiency, 
the semi-automatized process of preparing a code and the possibility of a calculation of 
the higher partial waves. The latter point gives hope for the future use of the full Tucson- 
Melbourne force in a description of 3N scattering at higher energies. The expected strong 
effects on observables coming from a 3NF should be tested also for short-range parts of 
three-body interactions. Such parts are included in the Tucson-Melbourne force. 

The automatized partial wave decomposition is especially important in view of future 
applications of 3NFs arising from the xPT. In this approach, consistent two- and three- 
body forces are derived The numerous spin- momentum and isospin structures, which 
occur at higher orders of the chiral expansion require an efficient and automatized method 
for the PWD. The here presented results for the TM force prove, that such a method already 
exists. 
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<p',q',a' I V^^\l+P) I p,q,a > [fm^] 
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FIG. 10: (Color online) The convergence of the full TM V^^\l+P) matrix elements with respect to 
the number of intermediate partial waves used during the action of the permutation operator. The 
channel combinations and momenta are as in Fig. [2j The crosses represent predictions obtained 
within aPWD approach. The dotted (black), dash-double dotted (red), dash-dotted (green) and 
solid (black) lines represent results obtained by the action of the permutation operator on the TM 
y(^) force with all the intermediate 3N states up to jmax = 2, 3,4 and 5, respectively. 
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FIG. 11: (Color online) The two-body correlation function for the triton for different NN potentials 
alone (thin blue lines) and together with the TM 3NF (thick black lines). Results obtained with 
the CDBonn, Nijmegen I, Nijmegen II and Nijmegen 93 potentials are represented by solid, dotted, 
dashed and dash-dotted curves, respectively. 

Appendix A: Standard PWD for it — p component. 

The ir-p part ^ of the Tuscon-Melbourne 3NF is given in terms of the momenta ki and 
ki ' of the individual nucleons as 



Introducing isospin, 1^ and /_, and spin, Fkr, Fj^j^, Fj^j^, Fj^_, and F^^_, operators we 
rewrite it as 



{k[k'^k'^\W^P\kik2h) 

= -^A^'i + k'^ + k'^-h-k2- h) 

xihRlW.l") FKR + I-RlMq\q") iFi,-F^,) + I^RZiq\q'') (F^. + F^.)} 



x{-(?fi ■ fa X f3)Rl^j^{q\q'^){iai ■ x ^ 
+ (r2 • mZi<l'^Q")i'l^ Q') ■ (^"'x ^2) 

+ ■ T2 X T:,)Rl'_{q\q'^)[{tai ■ x q')q^ - {la^ ■ qx q '){a2 ■ q)]} 
+ (2 ^ 3,^^ -q'). 




(Al) 



(A2) 



with 



(A3) 
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and 



f77 _ -{^3 ■ q') [iffi ■ q X g 0, 

qiz _|_ qz _|_ 



The R^x^i R^+, and form factors are given in terms of regularization form factors at 
the meson-baryon-baryon vertices Fi as [6|] 

RIW^ q") = ^i^p^^niq') + n,F,^NAq')]FpNNAq')F'NNiq") m 



and 



RZiq\q'") = \RZiq^q'") 

_ 1 gp m 5M — ^ 

{q')F^NA{q")F^NN{q"). (A9) 



The Fj are taken in monopole form 

FM') = (AlO) 

with z = {ttNNjTtNA, pNNi), pNNp, pNA}. The mass of the boson at the corresponding 
vertex, rrib, is either or rrip with exception of the case i = pNA when rrib = 0. 

We would hke to have matrix elements of the three-body force in a partial wave basis 
\pqa)i, where p and q are magnitudes of Jacobi momenta [p is the relative momentum 
between particles 2 and 3 and q is the momentum of spectator particle 1 relative to the 2-3 
pair) and a denotes discrete quantum numbers which we separate in spin, aj, and isospin, 
ax, parts 

\pqa)^ ^ \p{ls)jqiX^)lU, I)JM; {t^)TMT)i 
= Ws)3q{\\)I{jJ)JM),\aT)i 

= \pqaj)i\aT)i (All) 
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The partial wave states corresponding to different spectator nucleon i {i = 1,2,3) can be 
obtained from \pqa)i acting with proper permutation operator, for instance PisP23\pq<y)i = 
\p'q'a')3. 



According to the scheme presented in Ref. [16|] these matrix elements can be calculated 

poo poo 

,{pqa\I^Tnil + P)\p'q'a'), =V/ / p"'dp"q"'dq\{pqa\I^Tn\fq"a"h 

X3(pVV|(l + P)bV«')i (A12) 



with 



poo poo poo poo 

,{pqa\I±Tn\p"q"a")3 =V/ / p'"^ dp'" q'""" dq'" / p'^dpq'^dq 

^/// ^0 Jo ^ Jo Jo 

/•oo /»oo 

X / / p^dptfdq i{pqa\p"' q" a"')2 
^ Jo Jo 

X (2(/'/V;|J■(2)7^(2)|MdJ)2 2{maj\mM 
Xs{pqaj\J^^''^'R.^''^p'q'a'j)s) 2K\I±K)3, (A13) 

where spin operators J-" and form factors TZ are taken among F^r, F^+, F^+, F^_, F^L, 
and R^ji, and R^-, respectively, for different contributing terms. 



The matrix elements of the isospin parts appearing in Eq. (lASp are given [16| by 



2(aT|/-|ar)3 = 2((t^)™T|r2 ■ rs\{t'^)T' Mt')^ 



2 2 



STT'5MrM^,i-Q)i~Yyit'\ I 1 I (A14) 



t \ T 



A=l/2 k 2 2 ^ j ^/ 1 1 

' 2 2 



where we use abbreviation d = 2a + 1. 

In the following subsections we will present the resulting expressions for the partial wave 
decomposed matrix elements 2{p"'Q"'oi'j\J-'^'^^7l^'^^\pqaj)2 and s{P<io!j\J^^^^TZ^^^\p'(l'c('j)3 of dif- 
ferent contributing terms to tt — p Tucson-Melbourne 3NF. 



1. The KroU-Ruderman term FkrR'kr 

The matrix elements of 2{v"' q"' (^"J\J^^'^'^'R.^'^'^\pqaj)2 and 3{pqcij\J^^^^'R^^'^\p' q'oi'j)^ for the 
FkrR^kr term are identified as 

2{p'y"a"J\:F^'^n^^pqaj)2 ^ ^,2{p"'q"'a"J\ -^^'- ^V ,V(g^^)lM«.)2, 

lorn'' q'^ + 
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• — * — * ^ — * 

3(pgQ:j| ^2 1% [FpNNo{<f) + «p^piVJVp(g^)]i^pArAro(?^)|pVaj)3 
where q ' — p — p "' and q — p'—p- They are given by 

^P^^ / W /// ///I ~(^3 ■ ? ') 77.2 / /2n I . . . \ 
1^2(P g ^./l g/2+^2 ^^iViv(g )|P^«J)2 

_ ^^p^?' "... . 



and 



with 



X 



3 (PgO; J I ^^^^ ^ [FpNNn (g^) + l^pFpNNp{q'^)]FpNNn{(f)\p'Q'<^'j)3 

= ^/2 ^JJ'^AA'^/ 

x(-)W*'l2V67r\//''/' 



/' 1 / 
XI \V s's 





xlp<^-^^(p',p)-p'iff^-^^(p',p)), 





1 


I' s' j 






s i 1 




1 




(i 


1 


') 




2 






1 


s J 


2 


2 





PP \ ^^pNNp ^^pNNp, ^^pNNp ^^pNNn 



and 



-Dm. = T , -D — 



J^pNNp 



2p'p ' -^PiVAT^ 2p'p 
2p'p 



The are Legendre functions of the second kind. 
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2. The term i^A+^A+ 



The term is written as 



^+ «/2 _L m2 «2 I ™2 Z^>^ M V S «/2 I ™2 II «2 I ™2 r >^ ^ 



and the matrix elements of —fi component of 2(p"'g"'a'7|^*-^^^^^'*|pgQ<j)2 and fj, component 
of s{pqaj\J-'^^^TZ^^^\p'q'a'j)s for -^a+-^a+ identified as 



{2(/V'v;|^(2)7^(2)|pgd,)2} 

fp^ ^* mhM-m /47r g'Fi"^(g')(cT3 • g') „ . ,2n p. / /2m • • ■ \ 

is^^I^^-^MlW^"^^ ^ "-^'VT g-^ + m^ ^-""^^ ^^'^"^^^ 

3(pgaj| ^,^^^, [^pAriVg(g') + i^pFpNNp{q^)]FpNA{q^)\p'c^Oi'j)^, (A22) 



— * — * _ 

where q ' — p — p "' and q — p'—p- They are given by 

3 g'^ + 




5{q"' - q) 





1 

2 




J 










J'" 





2 1 1 

-407r^(-)""'+^' <i / s i 

s'" f" 



xJ2tHf^--^^{p"',p) Y: J^^^_ { r 1 ^i]c{aOlO,l"'0)CmO,iO) (A23) 



and 



g^ 

i{pqaj\al^— i^iFpNNoiq^) + i^pFpNNp{q^)]FpNA{q^)\p'q'a')3 

q^ + 



X I ], j, } C(1/.J'M', JM), (A24) 
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with 



p/l/p 



2 

m 



X 



'^'"'^ QKi?A..Jl, (A25) 



- 4- f^Ki^^J + -J^ QKi^A..J + ^-^^- 7- QKi?A..J 1 (A26) 



and 



p'p 



m 



2 



A2 



^ I f A2 _rr72VTO2-A2 1^'"*^^™"'* _ A2 VA^'' - A^ \Qi^^Knnb, 

+ 7T2 _ A2 VA2 T^^'^^'^^^p'va) 



(A^iVA--^)K-A^^^^)^'^^™^)+ (m2-A2^^^)(Af^^^-A2^J^^^^^^ 



^ -Qi{BApna) 



(ApATATp A^jy^)(A^^^ m 



(A27) 



The 5a^^^ and B^^na are given by 



n'"2 I ^2 I A 2 ^'"2 _|_ ^2 I A 2 

^A.^A - ^ , ^ApA^A - • (A28j 



The summation over fi in Eq. f lA2ip can be carried through resulting in 



- iiJM, JM)C{1 fiJ'M', JM) = 5jj'dMM'i~y~-^\l i. (A29) 
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3. The term 



The term is written as 



The matrix elements of — ;Lt component of 2{p"'Q"'ctj\J^^'^^T^^'^^\PQC(j)2 and component of 
s{pqaj\T^^^n^^^p'q'a'j)3 for are identified as 



AMr^^^^^Mir^""' ' "^"VT — — 



{sipqajlJ'^'^n^'^lp'q'a'j) 



3^ ^ 




47rgyi^(^)(a3-g) 



{PQ(^j\\l Y~L — 2 — ) + i^pFpNNp{,q )]FpNA{,q )\p'q a'j)3, (A31) 

q + ?7i 



where q' = p — p "' and q = p ' — p. The first term in Eq. ( ]A31I) is equal to Eq. ( 1A23I) and 
the second is given by 



AnqYi'mas-q) 




{pqaj\\ 5- 5 — [FpNNuiq ) + i^pFpNNp{q )]FpNA{q )bV«j)3 



3 q'^ + 771^ 
^S^_As^,.3^,.C(l ^J'M\ JM) 



r - - 

J 2 2 


1} 


{'r 


J 


5} 


\ 1 s' 


2 ■> 




J' 





X 



211 

-m7TV6{-y'+'^ lis] 



with 



I a+b=2 



H-np',p) 



V s' f 



(A32) 



( A2 ''^^2^KAnp) + J2 ^Z~V2 ^K^ApiViVo , 

PPY^pNA "V ^^pNNo ^^pNA 

A2 _^2 A2 _ A2 Vzl^ApivA. 
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A 2 a2 



PPy^pNA ^p ^^pNNp ^^pNA 

KnNp - KnA p 

^^pAfA ''V ^^pAfA ^^pNNp 



(A33) 



4. The term FI^R^^_ 



Using the identity 



/ 4-7r 

X a2 ■ g ' = zv^A/-g'5^(-)^a2^{ai,yi(g ')}''"^ (A34) 



3 



the term is written as 

{a^ ■ q ■ a2 ^ q ') q^ 



pi ^ 



r- ^^f {as-q')q'{ai,Y,{q')V'-^ U \ .A,,^ 

The matrix elements of — // component of 2{p"'q"'c(j'\J^^'^^T^^'^^\pqctj)2 and /i component of 
3(jigcij| J-''^^^??.*^'^^ |p'g'aj)3 for F^_R^_ are identified as 

i(/'g'"«';i^^'^^^'^|p^«j)2}"% 
1 9p9 m bM -m ^ 
4 48m5 *^^M M-m ""^ 

X2 P ^ "j v2W- ^- F^iVTV g )F^NA{q ) pqaj)2, 

V 3 g 2 -I- j7i^ 

)3} ^ 



3(Maj|^-7 — 2^FpNNniQ'^) + '^p^piVJVp(g^)]i^pJVA(g^)bVaj)3, (A36) 
g + ?Ti^ 

where g ' = p — p'" and q = p ' — p. The last term is identical to the last term in Eq. (1A22P 
for F^_,_ . The matrix element in the first term of Eq. ( ]A36p is given by 

2{P q aj V2\ — F^NNiq K^tva g ) pqaj)2 

V 3 g + 



yu ..111 



2 2 
2 2 * 
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a+6=2 



^ ^ li''' r } ^'"0)C(M/0, /O) 




(A37) 



5. The term FI^_R2'_ 



Using the identity 



a,-qx q' = -^^2^^^' ^(-j^^ja^, (A38) 



the term F^^_ is written as 

Fil ^ 



II _ -{^s ■ q') {^2 ■ q){i^i -qxq") 



= -V2Y2.(-r| ^,2 + ^2 II }• ^^39) 

J. II J. ^ 

The matrix elements of —jj component of 2{p"'q"'<^j'\J^^'^^'7^^'^^\pqC(j)2 and jj, component of 
3{pqdj\J^^^^7l^^^p'q'a'j)3 for Fj^_W^_ are identified as 



{2(p'"g'"a';|J•(^)7^(^)|p^d,)2}"' 
1 £/pf ^* mbM-m ,„ ,„ ,„ {^3 ■ q ')q'Y-''{q') 

-i^^^^^Mir^"^' ' "^'VY — ^^^T< — "^^^^^^ ^"^^^^^^ 

{3(pgci.|J-(=^)7e(=^)b'g'aS)3}'' 

V o ? ~r 

(A40) 

where q' = p — p'" and q = p ' — p. They are given by 
^^n{a3-q')q'yi'{q' 

S{q"' - q) 




2{p"'q"'a"J\\l^ '^' ' F^NN{q")F^r,^{q'')\pqajh 



3 q' +ml 



X 



X 



q 



u 


1 

2 














J'" 
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-407rv/6(- 



2 11 
I s j 
I'" s'" f" 



V(2a)!(26) 



(A41) 



and 



V O Q 'T' 



S{q' - q) , , 

^,2 ^A'A^/'/ 



X 



1 j s 



1 1 1 



1 1 

2 2 
lie' 
2 2"^ 



+2407rV6(-r^J](-)^fc|2 ^ J 

A; ^ 



2 A; 1 
/ s j 
I' s' f 



1 1 k 

1 i s 

2 2'* 
2 2"* 



X 



E 

a+b=2 



/a •■/) 



V(2a)!(26)! 



(A42) 



Note, that in the case of A^ttva = ^ttNN (as was assumed in [6| and also used by us in 
apphcations shown in this study), the function H[^'^~'^^^{p"' ,p) is equal to H[^^ {p'" ,p) 
defined in Eq. flATO]) and the function H'^^^-^^^{p'\p) is equal to Hf^^{p"',p) given by 



Hf'ip'^p) 



m 



p"'p 



72 {Ql{BmJ - Qi{BA^j^ff)) 



At 



NN 



m: 



- kl^^Q'j{B^^,,) (A43) 



2{p"'py 
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